Photoinduced Structural Phase Transitions in Polyacene 
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(Dated: 23 March 2011) 

There exist two types of structural instability in polyacene: double bonds in a cis pattern and 
those in a trans pattern. They are isoenergetic but spectroscopically distinct. We demonstrate 
optical characterization and manipulation of Peierls-distorted polyacene employing both correlated 
and uncorrelated Hamiltonians. We clarify the phase boundaries of the cis- and trans-distorted 
isomers, elucidate their optical-conductivity spectra, and then explore their photoresponses. There 
occurs a photoinduced transformation in the polyacene structure, but it is one-way switching: The 
trans configuration is well convertible into the cis one, whereas the reverse conversion is much less 
feasible. Even the weakest light irradiation can cause a transition of uncorrelated electrons, while 
correlated electrons have a transition threshold against light irradiation. 

PACS numbers: 78.30. Jw, 78.20.Bh, 78.20. Ci, 71.30.+h 



I. INTRODUCTION 

Since Salem and Longuet-Higgins pointed out that the 
Peierls instability in polyacene^ should be conditional 
and thus distinct from that in polyacetylene 2 -, its ground 
state has been studied by numerous theoreticians with a 
variety of numerical tools. Polyacene consists of linearly 
fused benzene rings and is the next simplest polymer to 
polyacetylene. Nevertheless, the order of the energetic 
stability in its structural isomers is still under discus- 
sion. The aromatic (AM) configuration of regular lat- 
tice closely competes with two types of Peierls-distorted 
structure: double bonds in a cis pattern and those in 
a trans pattern, which read as in-phase (IP) and out- 
of-phase (OP) bond order waves (BOWs) in a couple of 
rows, respectively. 

Standard- and extended-Huckel crystalline-orbital cal- 
culations at an early staged— pointed to the OP-BOW 
phase as the most stable. Band-structure calculations on 
a modified-neglect-of-diatomic-overlap (MNDO) basis 6,7 
followed to find the IP-BOW phase instead to be the 
most stable. While modern computational approaches 
such as quantum Monte Carle 8 and density-matrix renor- 
malization group^ also supported the IP-BOW phase 
favored over the OP-BOW one, there exist not a few 
arguments!^ that IP and OP BOWs are almost isoen- 
ergetic. They are indeed degenerate with each other un- 
less we take account of the electron itinerancy beyond the 
nearest-neighbor hopping^ and/or the electronic corre- 
lation beyond the Hartree-Fock (HF) schemed Various 
ah initio calculation !) further show that whether and 
how polyacene undergoes lattice distortion critically de- 
pend on the supposed electronic correlation and the im- 
posed boundary condition. There is thus every reason to 
believe that the cis- and trans-distorted isomers as well 
as the regular structure may coexist due to local defects 
and/or thermal excitations. 

Then a query, is there any idea of designing the elec- 
tronic structure of polyacene other than energetically? 



We demonstrate that optics can be the answer. We stim- 
ulate a renewed interest in polyacene: optical observation 
and manipulation of its isomers rather than energetically 
ranking them. 

There are transition-metal-based ladder materi- 
als and some of them2i~— may also be interest- 
ing in the context of varied structural instabili- 
ties. Among others are platinum-halide double- 
chain compounds, where charge-density-wave (CDW) 
states of the IP and OP types are both real- 
ized. (^-bpym)[Pt(en)Cl] 2 Cl(C10 4 )3-H 2 (/J-bpym = 
2, 2'-bipyrimidine = CsHgN^ en = ethylendiamine = 
C 2 H 8 N 2 )^ and (bpy)[Pt(dien)Br] 2 Br 4 • 2H 2 (bpy = 
4, 4'-bipyridyl = CioHgN 2 ; dien = dicthylcntriamine = 
C4Hi3N3)2£ were fabricated in an attempt to bring into 
interaction a couple of alternating platinum-halide lin- 
ear chains via bridging organic ligands, whose platinum 
ions are of mixed valence and form IP and OP CDWs, 
respectively i 26 ' 27 The intersite electron-phonon coupling 
brings BOW varieties into the hydrocarbon ladders, while 
the intrasite electron-phonon coupling brings CDW vari- 
eties into the platinum-halide ladders. Indeed the IP- 
CDW and OP-CDW states also exhibit distinct opti- 
cal features^ but they are never yet found in a sin- 
gle material. Because they are made in different en- 
ergy structures^ they are much less exchangeable with 
each other whether by light irradiation 2 ^ or by applying 
pressure^. On the other hand, the IP-BOW and OB- 
BOW states in polyacene are isoenergetic and thus able 
to coexist. Here is an increasing possibility of optically 
tuning optical properties. 



II. MODELING 

In order to elucidate the electronic correlation effect 
on the optical properties, we describe polyacene with 
and without Coulomb interactions, generally writing the 
Hamiltonian as 
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where cj ■ s and Ci: JjS (cj.^ s ci : j >8 = Hiy, s ) create and an- 
nihilate, respectively, a 7r electron of spin s =f,i= ± 
at site j on chain Z, while r; : j = — it; :J - signifies 

the bond distortion caused by the jth and (j + l)th 
carbon atoms on the Zth chain. We assign a stan- 
dard value, 2.4 eV^^ 9 - ' 31 ' 32 for the average intrachain hop- 
ping integral iy and work at a sufficiently low tempera- 
ture, feeT/tu = 0.01. If we set the Coulomb interac- 
tions all equal to zero, eq. ([T]) is reduced to the Su- 
Schrieffer-Heeger (SSH)-type Hamiltonian, 33 which is ef- 
ficient enough to illuminate the structural instability and 
excitation mechanism of polyacenej 12 ' 13 ' 34 " — In the SSH 
modeling, the electron-lattice coupling constant a and 
the c-bond clastic constant K are usually set equal to 
4.1eV/A and 15.5 eV/A 2 ^^^ respectively, together 
with the carbon-site mass M of 1350 eV • fs 2 /A 2 ^^ 
We further set the interchain hopping integral t±_ equal 
to 0.864t|| considering that C-C bonds in the rung (_L) 
and leg directions are made in different lengths as 
an = 1.4 A < a± = The full expression 

((T|) is regarded as an extended Peierls-Hubbard (EPH) 
Hamiltonian^ where we design the Coulomb param- 
eters as V||(j_) = U/n^Jl + 0.611701^42 k reads as 

a dielectric parameter and may range from 1 of bare 
correlatio n 31 ' 32 ' 43 to 2 of fully screened correlation . 32 ' 43 
Here we consider the intermediate case of K = 1.5 in an 
attempt to survey potential ground states and typical 
photodynamics in linear polybenzenoid hydrocarbons. 

We are interested in the generic behavior of polyacene 
rather than individual features of small oligoacenes. In 
view of the applicability of the periodic boundary condi- 
tion and the feasibility of simulating the electron-lattice 
dynamics on a large scale, we restrict Coulomb interac- 
tions to a certain range. We indeed learn from quan- 
tum Monte Carlo^ and density-matrix renormalization- 
group 9 calculations that distant Coulomb interactions 
have little effect on the ground-state properties, includ- 
ing the conditional Peierls instability and the energetics 
of the structural isomers. We have carried out tenta- 
tive calculations of both static and dynamic properties 
for 16, 32, 64, and 128 benzene rings. Optical features 



such as the absolute energy and relative intensity of every 
absorption are well converging even at N = 16. Quali- 
tatively the same dynamics is available as a function of 
the excitation density for N > 32. Therefore, we present 
typical findings at N — 64 unless otherwise noted. 

Figure Q] reveals possible ground states in polyacene. 
When the electron-lattice and electron-electron interac- 
tions are both weak, we find a metallic state of AM con- 
figuration. With growing coupling and correlation, the 
electrons are localized into an antiferromagnetic (AF) 
Mott insulator (MI) and a bond-alternating Peierls in- 
sulator, respectively. A CDW state also appears in the 
correlated region. Such a close competition between AF 
MI and CDW is due to the moderately valued screening 
effect. CDW is always favored over AF MI at n = l, 6 
whereas vice versa at k — 2^ Though the Peierls in- 
stability in polyacene is conditional, the well-established 
SSH parameters are known to result in a gapped ground 
state [cf. Fig. [2{a)]. The two types of Peierls insula- 
tor, IP and OP BOWs, are degenerate in energy. Even 
though Coulomb interactions are switched on, they re- 
main isoenergetic with their Peierls gap being enhanced 
[cf. Fig. HJb)]. Keeping the SSH parameters as they 
are, we turn the Coulomb interactions on. We discuss 
the correlated dynamics typically at U = 4.8 eV and in- 
quire further whether and how it changes getting toward 
a boundary to the CDW phase. 



III. SPECTROSCOPIC FEATURES OF THE 
ISOENERGETIC STRUCTURES 

Considering that a charge-transfer excitation costs en- 
ergy of eV order, we evaluate the real part of the optical 
conductivity as 



«r(w) = -^2\{Ei\J\Eo}\ 2 5{Ei-Eo-nu), (2) 



where \Ej) denotes the wave vector of the ith- lying state 
with energy E^ and J signifies the current operator, 
which is defined in the leg and rung directions, respec- 
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FIG. 1: (Color) A ground-state phase diagram on the a-U 
plain within the HF approximation. Crosses point to the SSH 
parameters without any correlation and the EPH parameters 
in the cases of U = 4.8 eV and U = 7.0 eV. 
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FIG. 2: Dispersion relations of the 7r-electron valence (A = 
1, 2) and conduction (A = 3, 4) bands based on the SSH mod- 
eling (a) and the U = 4.8 eV EPH modeling within the HF 
approximation (b). 
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3.0 4.0 

FIG. 3: (Color) Polarized optical conductivity spectra paral- 
lel (||) and perpendicular (_L) to the long axis for IP and OP 
BOWs, where every spectral line is Lorentzian broadened by 
a width of 0.09t||: Exact calculations based on the SSH mod- 
eling (a); HF (b), SCI (c), and SDCI (d) calculations based 
on the U = 4.8 eV EPH modeling. 
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The state vector is described in terms of Slater determinants and truncated at the order of double excitations: 
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|0) is the ground-state HF 



wave function with |0) being the true electron vacuum 
and creating an electron in the mth HF orbital with 
energy e m . The orbital counter m(/c,A, s) is composed 
of the wave number (k = 2nn/N; n = 1, • • • , N), band 
label (A = 1, • • • , 4) and spin quantum number (s = ±). 
Since an electron excitation is optically allowed without 
any momentum transfer, 44 the periodic boundary condi- 
tion serves to reduce the number of excited states to be 
considered. 

Excited states of the HF type consist of a single Slater 
determinant where f^°\i) = 8 0ti , f^\k, fj,,u, s;i) = 
Skfj,us,i and f {2) (k 1 ,k 2 , ^,^2,^1, <^2, si,s 2 ;i) = 0, giving 
Ei = hf(^'o|'H|£ , o)hf — e m{k,ii,s) + ^m{k,u,s)- Those of 
the CI type read as resonating Slater determinants, 47 
where their coefficients are determined so as to di- 
agonalize the original Hamiltonian (p}. Within 
the single-excitation CI (SCI) scheme, f^(i) and 
f( 2 \ki, &2, fii, fi2, v ii v 2, 8%, S2] i) remain unchanged and 
therefore only excited states are improved. Here is no cor- 
rection to the ground state, because no single-excitation 
Slater determinant is obtained by the Hamiltonian oper- 
ation on the HF ground state. Only configurations of the 
double-excitation type can directly interact with the HF 
ground state. Therefore, it pays most to truncate the se- 
ries expansion at the present order, which is referred to as 
the single-double-excitation CI (SDCI) scheme, though 
the cost of SDCI (~ N 6 ) is incomparably larger than 
those of HF (~ iV 3 ) and SCI (~ iV 4 ). 

Figure [3] shows comparative calculations of the opti- 
cal conductivity spectra for IP and OP BOWs, where 
only the SDCI calculation is carried out at N = 16. For 
the light polarized perpendicular to the conjugation di- 
rection, the two isoenergetic structures still look alike 
even in the SDCI scheme. The valence band of A = 1 
and the conduction band of A = 3 are both composed of 
molecular orbitals symmetric with respect to the reflec- 
tion about the plane bisecting every rung bond, whereas 
those of A = 2 and A = 4 consist of molecular orbitals 
antisymmetric to the reflection. Such properties of PM 
remain unchanged with bond distortion, whether the cis 
or trans type, while they do not hold in AF MI and CDW. 
Therefore, interband electronic excitations of the 1-to- 
3 and 2-to-4 types make no contribution to the rung- 
direction optical conductivity in both IP and OP BOWs. 
On the other hand, a long-axis-polarized photon clearly 
distinguishes between IP and OP BOWs. In the con- 



jugation direction, 2-to-3 interband excitations, includ- 
ing the highest-occupied-molecular-orbital (HOMO)-to- 
lowest-unoccupied-molecular-orbital (LUMO) transition, 
are optically forbidden on an IP-BOW background but 
allowed on an OP-BOW one. OP BOW can thus absorb 
long-axis-polarized photons of much lower energy than 
IP BOW can. The zone-center (k = tt) l-to-3 and 2- 
to-4 interband excitations are most contributive to the 
IP-BOW spectra in the conjugation direction, whether 
correlated or uncorrelated, while the OP-BOW spectra in 
the conjugation direction qualitatively vary with the rel- 
evant electronic correlation. Valence-to-conduction-band 
excitations of any type are optically allowed at the uncor- 
related level, but only near HOMO-to-LUMO transitions 
survive the strong correlation. 

IP and OP BOWs may be coexistent in actual 
oligoacenes. They both exhibit selective absorption of 
their own for the light polarized parallel to the conjuga- 
tion direction, especially under the influence of the elec- 
tronic correlation. Then it is highly interesting whether 
and how the isoenergetic structural isomers are photoac- 
tive. Sufficiently oriented long acenes may have an ap- 
plication to the optical memory. 



IV. PHOTOINDUCED PHASE TRANSITIONS 
BETWEEN THE ISOENERGETIC STRUCTURES 

We simulate time evolution of a charge-transfer exci- 
tation solving the Schrodinger equation 



im s (t)=n s (t)v s (t), 



(6) 



where the spin-s sector of the Hamiltonian, T-L s {t) 1 and 
a complete set of its wave functions, & s {t), are given by 
square matrices of degree AN. The SSH Hamiltonian 
is exactly computable, while the EPH Hamiltonian is 
treated within the HF approximation. In both cases, 
defining the time-dependent wave functions as 
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eq. © is explicitly written as 
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where the electron-density and bond-order parameters 
are calculated through 

di-.j(t) = EE' l^w(*)| 2 , (9) 

S Si 

plj(t) = E/'/,,- ,:/i =EE'4,.( t )fe+w(i). (io) 

S St 

s s z. 

with denoting a summation over the occupied levels 
immediately after the photoexcitation. 

Discretizing the time variable as t m = mAt (m = 
0, 1, 2, • • •) with an interval much smaller than the optical- 
phonon time scale, say, At = 10~ 3 -J M j K ~ 0.01 fs, we 
integrate eq. ([6]) step by step: 



& s (t m +x) = cxp 



-At 



& S {t m ). (12) 



At t = to, the system is photoexcited, that is, some va- 
lence electrons are pumped up into a conduction band 
across the Peierls gap. In response to the photoinduced 
interband electronic transition, we immediately construct 
the initial Hamiltonian H s (to) from the modified elec- 
tronic parameters di-j(t a ), pj. (to), and pf{t a ) leaving the 

background lattice untouched. Then lf - s (to)^'H s (to)^'s(io) 
is no more diagonal and that is why the wave vectors be- 
gin to fluctuate. Obtaining the instantaneous eigenfunc- 
tions and eigenvalues 
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from H s (t m ) at every time step,— we can express eq. 
©as 



47V 2 2N 
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(15) 



Once the electronic wave functions deviate from the 
ground-state equilibrium, the lattice begins to fluctuate 
correspondingly. Its dynamics is governed by Newton's 
equation of motion, 

Mill:j(t) = Klm-.j-xit) - 2ui :j (t) + u i:j -i(t)] 

+2aRe[pl(t) - pj^^t)] = F l:j (t), (16) 

and solvable stepwise: 

Ul:j(t m+ l) = UUj(t m ) H — /It, 

W{:j(*m+l) = Ui:j(tm) + Ul,j{t m )At. (17) 

The classical treatment of lattice degrees of freedom 
is successful in investigating both correlated^ and 
uncorrelated^S electron-lattice dynamics. We avoid set- 
ting any "artificial" fluctuation on the lattice initially but 
start every calculation from a stationary BOW lattice of 
the IP [a = 1) or OP (a = -1) type: 



<J l ri.j{t ) = -a l ri :j+1 (t ), ui : j(t ) = 0, 



(18) 



in an attempt to simulate every precursor to structural 
phase transitions quantitatively. 
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FIG. 4: (Color) SSH-model calculations of the bond variables r± : j and r; :j as functions of space j and time t in the case of the 
HOMO-to-LUMO one-electron excitation on the OP-BOW background. 



120 - 
100- 
B0- 

;' ■ 

60- 
40- 
20- 



I* 



0.0 



120- 
100;- 

30 

J 
60 

40 

20 



1.0 f(ps) 2.0 
'I'' 



MM 




120 
100 
80 



SO 
40 



r 2:J 



-0.5 



0.0 



0.5 



1.0 t(pn) 2.0 



3.0 



■ J I 



120 
100 
SO 
-I 60 
40 
20 



FIG. 5: (Color) The same as Fig. [4] but in the case of a 2-to-3 interband three-electron excitation at k ~ n. 



In order to visualize bond arrangements of the IP 
and OP types, we introduce effective bond variables, 
(ji-.j ± r2y)/2 = r±;j, where the bare bond distortions 
r\.j are decomposed into their net (fi : j) and alternating 
(fUj) components as fi-j = (2>/ : j + ri-j-i + r/ :j *+i)/4 and 
ri-.j = (2ri : j — ri-.j-i — r /:j+i)/4. ri-.j remaining constant 
with varying j read as regular bond alternation on the ith 
chain, while a change of sign in them denotes a kink (do- 



main wall) occurring. 



r+:j 



and 



serve as detectors of 



the IP-BOW and OP-BOW configurations, respectively. 
Now we are ready to observe photoexcitations of poly- 
acene. We consider, unless otherwise noted, a few elec- 
trons being excited by around Peierls-gap energy, that is, 
transitions from the highest one or more occupied to the 
lowest one or more unoccupied molecular orbitals. The 
thus-induced electron-lattice dynamics depends on the 
momentum and spin of the excited electrons but exhibits 
universal features. 
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FIG. 6: (Color) The same as Fig. gjbut on the IP-BOW background. 
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FIG. 7: (Color) The same as Fig. [5] but on the IP-BOW background. 



A. Uncorrelated Electron-Lattice Dynamics 



We begin with correlation-free photoresponses. In 
Fig. @]we demonstrate the most fundamental HOMO-to- 
LUMO one-electron excitation of the OP-BOW ground 
state, which is optically allowed in both leg and rung 
directions. The photon energy is first spent on a back- 
ground uniform oscillation and then contributes to an 
OP-to-IP-BOW phase transition. The two bond-ordered 
chains are disarranged at the same time. The defects in 



both chains split into a pair of kink and antikink. Then 
the resultant antidomains grow "repulsively" . The out- 
of-phase lattice fluctuations convert the OP-BOW ar- 
rangement into the IP-BOW one. The OP-BOW con- 
figuration macroscopically recovers within a picosecond 
but declines again, and finally disappears into the global 
IP-BOW configuration. Stronger light irradiation drives 
more electrons to jump over the Peierls gap, induces 
a precursory global oscillation of larger amplitude and 
longer period, and then generates multipairs of kink and 
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antikink against the initial bond arrangement. Figure 
[5] shows the lattice relaxation from a low-lying three- 
electron excited state. Declining and recovering again 
and again in such a complicated manner, OP BOW dis- 
appears into IP BOW within a few picoseconds. Such 
a swinging motion is observed in a double-chain CDW 
system without any Coulomb interaction as well^i and 
therefore recognized to be characteristic of uncorrelated 
electrons. 

What happens to IP BOW being photoexcited in the 
same way? Figure |B] demonstrates the HOMO-to-LUMO 
one-electron excitation of the IP-BOW ground state and 
the following lattice relaxation. Indeed the absorbed pho- 
ton similarly brings a pair of kink and antikink to both 
chains, but the resultant antidomains now grow "attrac- 
tively" . The kink-antikink breathing motions in the two 
chains are completely in phase and there is no sign of 
OP BOW appearing. IP BOW withstands stronger light 
irradiation. Figure [7] shows that IP BOW with three elec- 
trons being excited still keeps itself within a configuration 
of the IP type. Pumping further electrons up to the lower 
conduction band, we find no sign of OP-BOW nucleation, 
even though numbers of antidomains eat away into the 
initial configuration. Thus excited IP BOW is expected 
to relax, whether radiatively or nonradiatively, into the 
initial state via kink-antikink geminate recombinations. 

We have so far investigated low-energy excitations of 
near Peierls gap. OP BOW can be converted to IP 
BOW by the weakest light irradiation, while IP BOW 
is light-resistant and survives much stronger photoirra- 
diation. However, there is a possibility of such an in- 
active response of the IP-BOW state resulting from the 
light polarization. We should be reminded that a long- 
axis-polarized photon is not contributive to any inter- 



band electronic excitation of the 2-to-3 type on an IP- 
BOW background (cf. Fig. E). Therefore, Figs. [Hand 
just say that charge-transfer excitations in the rung di- 
rection hardly lead to an IP-to-OP bond-structure con- 
version. We should inquire further into a possibility of 
phototuning IP BOW by higher-energy excitations. We 
thus try the 2-to-4 interband one-electron excitation at 
k = 7r, which reads as the optically-allowed lowest-energy 
charge-transfer excitation in the conjugation direction on 
an IP-BOW background, but the resultant lattice mo- 
tion is almost the same as Fig. [5] without any sign of 
OP-BOW configuration appearing. When we excite fur- 
ther electrons to the upper conduction band (cf. Fig. 
[5]) , two chains are no more exactly harmonized with each 
other. Transient islands of OP-BOW configuration come 
to appear here and there with increasing number of elec- 
trons being excited, but they do not grow up completely. 
IP BOW looks so persistent and tough against any light 
irradiation. 



B. Correlated Electron-Lattice Dynamics 

Now we switch on Coulomb interactions and discuss 
correlated-electron-driven lattice dynamics in compari- 
son with the above. When we excite one electron across 
the Peierls gap on the OP-BOW background in the same 
manner that we did in Fig. 01 the ground state remains 
silent without any explicit domain-wall nucleation (Fig. 

With increasing number of absorbed photons, OP 
BOW happens to change into IP BOW with three elec- 
trons being excited (Fig. ITU)) but fails to do so even with 
four electrons being excited (Fig. [TT1) . OP BOW seems 
to be necessarily converted into IP BOW with more than 
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FIG. 9: (Color) U — 4.8 eV-EPH-model calculations of the bond variables r±-j and ri-j as functions of space j and time t in 
the case of the HOMO-to-LUMO one-electron excitation on the OP-BOW background. 
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FIG. 10: (Color) The same as Fig. [9] but in the case of a 2-to-3 interband three-electron excitation at k ~ n. 



four electrons being excited. When electrons are corre- 
lated, an OP-to-IP conversion is not feasible until the 
light irradiation comes up to a certain threshold intensity. 
Even when an excitation density is given, the correlated 
dynamics still depends crucially on how electrons are ex- 
cited. The initial excitation density is determined by the 
light intensity and a threshold in it, if any, should origi- 
nate from an energy barrier in the relaxation path. The 
critical number of excited electrons is not so decisive as it 



is in a single chain^ because the spatial degrees of free- 
dom in charge-transfer excitations due to the ladder-like 
lattice structure of polyacene may bypass the relaxation. 
That is why the threshold excitation density is signifi- 
cantly influenced by the microscopic information about 
photocarriers such as their momenta and spins. 

What about photoexciting IP BOW under the influ- 
ence of Coulomb interactions? Figure Q2] shows the case 
of three electrons being excited. The observations re- 
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FIG. 11: (Color) The same as Fig. [9] but in the case of a 2-to-3 interband four-electron excitation at k ~ tt. 
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FIG. 12: (Color) The same as Fig. flOlbut on the IP-BOW background. 



main qualitatively unchanged from the uncorrelated case. 
It indeed happens that local OP-BOW domains appear 
with further electrons being excited. However, they 
hardly aggregate to develop a long-ranged order and de- 
cline within several picoseconds. Though the electronic 
correlation may stimulate local OP-BOW domains to ap- 
pear, yet a global IP-to-OP-BOW phase transition still 
looks hard to photoinduce. 

Finally we increase the Coulomb interactions so as 
to get closer to a phase boundary to CDW and pho- 



toexcite IP and OP BOWs by around Peierls-gap en- 
ergy. Figure Q2] (THl) demonstrates the electron-lattice 
dynamics with four electrons being excited from the OP- 
BOW (IP-BOW) ground state, where we define an ef- 



fective site variable, di : 



d 2 : 3 = d- 



■j • 



with duj — 



(2di : j — di.j-i — e£; : j+i)/4 and plot it as well as r± : j in 
order to detect the CDW configuration. In both cases, 
electrons ultrafast form up into the CDW order and then 
the BOW lattice arrangement melts down. The liquidized 
BOW state should fade out if we include the effect of 
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FIG. 13: (Color) U = 7.0 eV-EPH-model calculations of the 

bond variables r± : j and the site variable d-j as functions 
of space j and time t in the case of a 2-to-3 interband four- 
electron excitation at k ~ 7r on the OP-BOW background. 
The initial 50 femtoseconds are magnified in an attempt to 
visualize the ultrafast dynamics in more detail. 



extrinsic dissipation into our calculation. A global and 
surviving CDW order is not attained until the photoex- 
citation density reaches a certain threshold. The devel- 
opment of the CDW order simply results from closing 
in upon the phase boundary rather than increasing exci- 
tation energy, because such observations are hardly ob- 
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FIG. 14: (Color) The same as Fig. [13] but on the IP-BOW 
background. 



taincd from every other previous parametrization. We 
have indeed tried various photoexcitations, including in- 
terband electronic transitions of the 2-to-4 type, in both 
SSH and U = 4.8 eV EPH models, but such a global 
CDW formation is of no occurrence. It is not the BOW- 
to-CDW transition itself but lattice lagging behind elec- 
trons in their combined dynamics that interests us here. 
In this context we may be reminded of a charge-transfer 
compound, tetrathiafulvalene-p-chloranil, where the elec- 
tronic valence instability and the dynamic lattice insta- 
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bility exhibit different time scales <£2 First an ionic do- 
main is photogenerated, whose rising time is so fast as 
to be less than the experimental resolution, and then the 
constituent molecules are stimulated to oscillate in a spin- 
Peierls-like manner. Both observations can be charac- 
terized as charge transfers preceding lattice oscillations. 
There is particularity in addition in the present findings. 
Even in such ultrafast dynamics IP BOW gives an indi- 
cation of its superiority over OP BOW. OP BOW is once 
rearranged into IP BOW and then liquidized, while IP 
BOW directly breaks down without explicit nucleation 
of OP BOW. 



V. SUMMARY AND DISCUSSION 

We have optically characterized and manipulated the 
two isoenergetic structures of polyacene with bond ar- 
rangements of the IP and OP types. An OP-to-IP-BOW 
phase transition can be induced by the weakest light ir- 
radiation unless electrons are correlated. Coulomb in- 
teractions function as energy barriers in the relaxation 
path and then the transition is conditioned on a good 
number of electrons being excited. On the other hand, 
an IP-to-OP-BOW phase transition is much less caused 
by photoexcitations, whether electrons are correlated or 
not. IP and OP BOWs may coexist in polyacene. If 
the Peierls-distorted polyacene is irradiated with light of 
around Peierls-gap energy, the mixed ground state is pos- 
sibly unified into a monostructure of the IP-BOW type. 
The photoinduced predominant state can be annealed 
back into the bimodal configuration. Figure [3] shows 
that IP and OP BOWs selectively absorb ultraviolet and 
blue-green rays, respectively, provided ti\ — 2.4 eV. Poly- 
acene thus reveals itself as a photochromic polymer 53 
with novel structural bistability. 

Photoinduced charge-transfer excitations potentially 
develop into kink-antikink bound pairs. Their breath- 
ing motions can cause a global phase change. Figure [TTJ 
demonstrating geminate recombinations of domain walls 
in the early stage, implies that correlated electrons en- 
counter energy barriers on the way of photogenerated 
excitons being self-trapped and further dissociated into 
distant kinks and antikinks. A correlation-induced en- 
ergy barrier on the relaxation path to far distant domain 
walls is observed in platinum-halide chain compounds as 
well ) 46 ' 48 ' 54 ' 55 which are well known to be photoactive 
Peierls insulators. Since a valley precedes the barrier, 
Coulomb interactions stimulate the Frank-Condon state 
to relax and accelerate its localization.™ We are thus con- 
vinced of the ultrafast quasifemtosecond self-localization 
of photogenerated charge-transfer excitons in the corre- 
lated dynamics (Figs. ITUHT41 and its slowing down in the 
uncorrelated dynamics (Figs. HHHJ)- 

Coulomb interactions further serve quick and certain 
switching. Comparing Figs. [S] and [TU1 we learn that 
the electronic correlation suppresses the macroscopic os- 
cillation of going back and forth between IP and OP 



BOWs. In order to illuminate such a contrast between 
correlated and uncorrelated electrons in more detail, we 
try to quantitatively evaluate the weights of IP-BOW 
and OP-BOW configurations. Their time-dependent and 
spatially- varying amplitudes may both be normalized in 
comparison with the initial uniform BOW, because co- 
existent stable IP-BOW and OP-BOW solutions have 
exactly the same amplitude. In the case of photoexcit- 
ing OP BOW, the normalized momentary strengths of 
IP-BOW and OP-BOW configurations can be defined as 

P±(t) = Ef=i P±..j (*)/ Eli rL-.M We calculate P±(i) 
with and without Coulomb interactions and plot them in 
Fig. US] Besides the precursory global oscillation without 
any OP-to-IP-BOW conversion yet, the initial OP-BOW 
state composed of uncorrelated electrons repeats rise and 
fall on a macroscopic scale for a picosecond or less. The 
uncorrelated OP-BOW state is so persistent in its config- 
uration. On the other hand, the initial OP-BOW state of 
correlated electrons is much more quickly and smoothly 
replaced by the photoinduced IP-BOW phase once the 
excitation density reaches the threshold. The promotion 
of OP-to-IP-BOW conversion is also recognizable as a 
correlation effect. With a few electrons being excited, the 
converted fraction amounts to four fifths in the correlated 
case but stays less than two fifths in the uncorrelated 
case. Correlated electrons possess an adiabatic potential 
energy surface full of ups and downs. Once they climb 
over the energy barriers and find a stable configuration, 
they are unlikely to swing back to the initial state. 

The macroscopic oscillation of uncorrelated electrons 
in their photoconversion to the IP-BOW state and the 
correlation-driven globalization of the IP-BOW domain 
are both understandable by observing kink excitations 
in more detail. The phase transition originates in gen- 
eration of kink and antikink in pairs. Within the SSH 
modeling, neither spin nor charge accompanies any kink 
as a result of resonating nonlinear excitationS ) 48 i 56 where 
neutral and charged kinks are isoenergetic. With increas- 
ing Coulomb interactions, their degeneracy is lifted and 
those with excess spin (charge) density in their centers 
are relatively stabilized (destabilized). In the EPH cal- 
culations some of the kink excitations indeed bear net 
spins and it is their intrachain interactions that develop 
and hold the induced IP-BOW arrangement against the 
OP-BOW background. In a successful OP-to-IP transi- 
tion [cf. Figs. El UM C '), andQSKd')], we find two or 
more kinks conveying net spin densities of the same sign 
in each chain, whereas in a frustrated one [cf. Figs. |H1 
ITST a'). and [TSTb'')]. we observe bound pairs of kink and 
antikink with excess charge, rather than spin, densities 
disappear in the early stage. The repulsive same-spin- 
carrying kinks may be the key to a global phase change. 
On the other hand, SSH kinks without any excess density 
of spin and charge look hardly correlated and can there- 
fore get closer to each other (cf. Fig. |4|). They are more 
and more created with increasing number of photons be- 
ing absorbed, but their geminate recombination is of fre- 
quent occurrence. That is why the converted fraction P + 
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FIG. 15: (Color) Probabilities of finding the IP-BOW and OP-BOW configurations, P±, as functions of time t after photoexci- 
tation: SSH-model [(a) to (d)] and U = 4.8 eV-EPH-model [(a') to (d')] calculations in the cases of 2-to-3 interband one-electron 
[(a) and (a')], two-electron [(b) and (b')], three-electron [(c) and (c')] and five-electron [(d) and (d')] excitations at k ~ ir on the 
OP-BOW background, (a) and (c) correspond to Figs. [4]and[5] respectively, while (a') and (c') to Figs. [9land [l0l respectively. 



oscillates on a macroscopic scale in the SSH modeling [cf. 
Figs. [FoTa)- [T5l c) in particular]. The reduction of both 
IP-BOW and OP-BOW strengths with increasing excita- 
tion density [Figs. [r5la) 4T5T d)] is also attributable to the 
free kinks and thus-unlimited domain-wall nucleation. 

A variety of photoinduced phase transitions have so 
far been observed and analyzed indeed. However, most 
of them assume the initial condition to be close to the 
phase boundary and such a situation often depends in 
quite a critical manner on some structural parameters 
which are hardly tractable. Oligoacenes possess inher- 
ent bistability of topological origin and their Peierls- 
distorted structural isomers are strictly and stably isoen- 
ergetic. Even under uniaxial pressure, for instance, ap- 
plied in an attempt to tune their intrinsic properties such 
as t±/t\\ and Vj_/Vji, the strict bistability holds and en- 



ables us to systematically investigate whether and how 
energetically degenerate isomers are photoconverted into 
each other. Polyacene derivatives and analogues such as 
polyphenanthrene^ and polyacenacene^. will make possi- 
ble further comparative studies. We are looking forward 
to more and more interest in this potential organic poly- 
mer. 
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